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The general characteristics and the bandwidth response 


of the Kaijo Denki PAT 313-1 sonic anemometer/t 


hermometer 


were examined. Measurements were compared to those of a 


hot wire anemometer system and a one mil thermocouple sys- 


tem. Turbulence data were obtained by all three systems 


simultaneously in two experiments conducted at 
site. The sonic anemometer/thermometer system 
be a reliable, rugged, easily installed system 
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I: TINERODUCT EON 


A: PROBLEM 

Turbulent processes in the atmosphere are of interest 
to meteorologists since it is by these processes that effi- 
cient vertical and horizontal transfers of energy, heat 
and momentum occur. These transfer processes have been 
experimentally observed to occur within the frequency band 
OOO AZ to 10.0 Ha. An inportant method for describ- 
ing the role of turbulence is tnrough use of a system of 
instruments capable of measuring over this band. The ana- 
log data obtained is then read into computer systems to be 
digitized into discrete samples and then to be analyzed by 


spectral or other techniques. 


Pee OBJECTIVES 

The primary objective of this thesis was to utilize a 
sonic anemometer/thermometer to measure SE 
evaluate the band over which it is capable of measuring. 
A secondary objective was to utilize a Hybrid computer sys- 
temto digitize the analog data and then to use the IBM- 
360/67 computer system for spectral analysis of the digi- 
tized data. A third objective was to compare a Do 
anemometer/thermometer with a hot wire anemometer system 


and a one mil thermocouple system. 
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PI, GENERAL DESCRIPTION AND OPERATING PROCEDURES 
OF THE 
SONIC ANEMOMETER/THERMOMETER MODEL PAT-313-1 
A. GENERAL 
Tne Ek Bene nermometer- measures the three 
velocity components of the wind by detecting accurately 
the transit time difference (transit time sum for tempera- 
ture) between two acoustic pulses propagated in opposite 
directions across a fixed path. The three velocity com- 
ponents of the wind are the longitudinal, lateral and ver- 
tical components denoted by u, v, and w respectively along 
axes X. Y and Z (Figure 1). The w component is obtained 
directly from the vertical axis, but u and v have to be 
computed from the velocity components U and Va measured 


along the horizontal sonic paths A and B which are 120 


degrees apart. 


BE COMPONENT PARTS OF THE SONIC ANEMOMETER/THERMOMETER 
MP ormer Supply Indicator Unit (PR-1) 

The power supply/indicator unit and its component 
breakdown is pictured in Figures 2, 3, 4. The instrument 
operates on 100/115 volts with + 10% variation at 50 to 
60 Hz. Total power required is 100 watts. All switches 
and controls are located in this unit. Calibration and 
adjustment procedures (discussed later) are performed at 


oe PAT-313-1 manufactured by the Kaijo Denki Company 
Eimrtegsof Tokyo, Japan. 
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a Min meter output (analog form) is available 


era ndemeasurement on three range scales of 1) O to + 1 


ieee, 2990 vom Som/sec and 3) 0 to + 10 m/sec (Fig- 
ure 3). Temperature meter output is + 10C from mean air 
temperature which is dialed in on 9 dial (Figure 4). Two 


BNC output jacks are available for tape recording or sig- 
nal monitoring purposes on both wind and temperature units. 
2. Probe (TR-31) 

The probe has to be employed within the wind direc- 
mI IImits indicated in Figure 5. The wind velocity would 
be less than actual if the probe were positioned so that 
the wind was parallel to the path travelled by the sound 
pulses or if the supporting structure of the probe inter- 
rupted the wind flow. 

The probe heads are ceramic type transducers from 
which sound waves are trasmitted at a frequency of 100KHz. 
The transducers have high directivity and damping to mini- 
mize interference with other channels and are small physi- 
pM order to avoid blocking the air flow (Figure 6). 
The transducers screw into special holders designed to 
minimize acoustic coupling through the array. All probe 
heads are interchangeable with only minor adjustments re- 
guired in the wind/temperature units affected by the 
change. 

The probe heads, holders and entire array structure 
have been designed to reduce interference to the wind flow. 


All wiring and connections have been inserted inside the 
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A,B— axes along sonic paths. 


X,Y— axes of symmetry for 
: the arrays 


X E J 
sy horizontal vector=mean 
axes ol the- wind. 


USEFUL WIND 
DIRECTION 
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ripe 1. Wind and Sonic Axes. 
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POSI for protection from the elements as well as 
e educ inge interference (Figure 7). 
3. Junction Box (0A-33H) 

The junction box contains three sets of pulse os- 
cillators (which generate the LOOKHz signal sent to the 
transducers) and six sets of head amplifiers (which amplify 
by 1000 the 100KHz signal sent to the transducers). The 
wnet ron box is a watertight unit with connectors at both 
ends (Figure 8). One hundred meters of cable link the 
junction box and supply/indicator unit and three meters of 
cable link the junction box and probe (Figure 9). The 
ee ton bor must be connected to the probe and power sup- 
ply prior to turning the sonic anemometer on or the output 
transformer will be shorted. 

C. ADVANTAGES/DISADVANTAGES OF THE SONIC ANEMOMETER/ 

THERMOMETER 

The sonic anemometer/thermometer has many advantages 
over previously developed wind and temperature measuring 
systems. Recent PEE by McBean (1972) revealed the 
advantages the sonic anemometer has over a Gill propeller 
anemometer especially with respect to sensor response and 
angle of .attack (cosine response) in eddy correlation 
measurements. A report by Kaimal et al (1971), in which 
they compared the pulse and continuous-wave type sonic 
anemometers, lists many advantages of the sonic anemo- 


meter system over other velocity measuring systems. 
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cb) Side view 


(a) Top view 


A-comrcnent 
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Wind Direction Limits of Probe Type TR-31, 
Kaijo Denki (1969) 


Figure 5. 





Transducer 
( 79 /00kKHZ Ceramic) 


Figure 6.  Probe Head.  Kaijo Denki (1969) 
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Breure 7. 


20 








Figure 9, 
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The PAT-313-1 system has no mechanical parts which 
eliminates lag due to mechanical inertia and thus allows 
the sonic to measure very light and sudden changes of the 
wind. The sonic also has high sensitivity and a broad 
range in all three wind components and in the temperature 
component. 

The resolution of the sonic anemometer is t 2 cm/sec 
in a span of 0-10 m/sec and the resolution of the sonic 
Bu srspomerter is + O0.1C in a span of + LOC centered around 
a mean temperature range of -10C to 40C. The response is 
linear and so is its sensitivity. (Therefore, the res- 
ponse is the same at 5 cm/sec as it is at 10 cm/sec.) 

The sonic anemometer has a time response in the milli- 
second range as it detects changes in pulse time between 
two sets of transducers. It takes measurements 440 times 
[sec for the wind components and 220 times/sec for the 
temperature component which allows it to detect changes of 
approximately 100 Hz at an accuracy of three percent 
(Kiajo Denki, 1969). 

The sonic anemometer is not affected by exposure to 
the direct rays of the sun or to the rain as would a thermo- 
meter or a hot wire anemometer. It is, however, affected 
by immersion in water. A high moisture, high heat environ- 
ment affects measurement procedures. These have to be 
corrected in order to get true readings because the measur- 
ing circuits are set for an atmosphere of zero percent 


relative humidity at a temperature of 20C. The probe, 


ER 





junction box and connecting cables are waterproofed allow- 
ing for all weather operation of this system. There are 
lotes loe sensitive elements (such as a hot wire) 
exposed to the weather. It is a rugged system that can 
accept some rough handling. Cable connection points, es- 
pecially where they enter the probe, are sensitive areas. 
Too much twisting of the cables can result in broken wires 
in the connectors. Care must be exercised so as not to 
bend the probe thus changing the length of transmission of 
the acoustic path and causing misalignment of the probe 
heads. Both problems would create serious errors in velo- 
city and temperature een NN Each of the probe heads 
has a rubber cap over it to protect the ceramic transducer 
element but abrasive handling will rip it. 

The probe has a small bubble level in the top plate 
which is used to insure that it is properly leveled when 
installed. Errors are introduced if the probe is not hori- 
M The probe should be mounted in such a way that 
it is free of vibrations. It measures relative speed of 
the environment with respect to the probe so any movement 
of the probe will introduce errors into the signals. Any 
object placed within two meters of the probe will reflect 
erroneous signals into the transducers. These are recog- 


nizable as noise spikes. 





Kaimal et al (1971) has an excellent discussion of these 
errors, 


ge: 
Anat a 1971) discusses this error. 
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PERI ASUREMENT PROCESS 
1. Wind Measurement 

As described previously, the sonic anemometer mea- 
sures wind velocity by detecting the transit time differ- 
ence between two acoustic pulses propagated in opposite 
directions across a fixed path. The nomenclature used to 
describe the measuring process for a typical pair of 
transducers is shown in Figure 10. 

The propagation times (ti, to) for the pulses from 


TR, to Ri and TR, to R, are 


si u C cos a + V (1) 


and 


L 
m GC cos a — Va (2) 


The travel time (t4) of the sound with the wind is less 


Aman the travel time (t,) with the sound against the wind. 


The difference At = to 7 €y is given by 
ZEN, 
EE (3) 
(C cos a) - ES 


Angle a is assumed to be small; therefore, cos a 





% 1. One further assumes that Da > v? and therefore, 
Z 
C I> > De These approximations result in 
A (4) 
2 d 
C 
which, when solved for Va gives 
2 
NEG CAE 
EE or (3) 
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V Wind Velocity 
and 
V Sound Transmitting 


K 
KEY: » distance between transmitter and receiver 
C sound velocity 
V wind velocity 
VE component of wind velocity perpendicular 
to sound transmitting path 
Va component of wind velocity parallel to 
sound transmitting path 
t mentor sound to travel from YT, to R 
1 A e 3 Bt 
(with wind) 
t Cimemuome SOMmMdNEO.Eravel from T, to R 
2 E 2 2 
(against wind) 
IR), IR, transmitting transducers 
Ri, R, receiving transducers 
Q angle between transmitting path ana dis- 
tressed sound path 
Figure l0.  Schematic of Typical Measuring Process and 


Nomenclature. 
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The wind velocity parallel to the sound transmitting path 
(v) is therefore determined and is the guantity indicated 
on the metered output of the sonic anemometer as U, Or Ve: 
Nieva lues ot au and v (the longitudinal and lateral wind 
components of the X and Y axes of Figure 1) have to be com- 
puted from the velocity components U, and VI The w com- 


ponent is taken directly from the output. 


The equations for u and v are: 


u= (cos ß - i sin’ 6) U, + (cos ß + ES sin f) Va (6) 
3 v3 

1 - 1 å 
v =-( — cos B + sin BR) UT ( — cos B - sin E (87) 

Y3 /3 
where " Nu 

V - U 
ze tan”? B A 
/3 GA E U,) 


Overbar denotes time average over a suitably long 
period of time  Kaimal et al (1971), 
ES Temperature Measurement 
Temperature is measured in a manner similar to the 
wind measurement except that the propagation times (t, t5) 
are added instead of subtracted. The w component trans- 
ducers provide the data used for this measurement. The de- 


fining expression for temperature is 


24 (c? E UN EIE 
C - e 
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Again assuming c? NS vi and c? ps v? 
eguation (8) reduces to: 
2% 
SE co) 
The speed of sound in air is expressed as: 
AO e +0.32 £) 17? (10) 


where T is in degrees absolute (K), e is the vapor pressure 
in atmospheres and p is the atmospheric pressure in atmos- 


pheres. From (9) and (10) 


20027 
Sore 
Je E EM OU a (LI) 
zi E : 
Furthermore, set Toy = T(1 + 0.32 SS where Toy is 


sound virtual temperature (taken as the true temperature 
in degrees K) measured by the sonic thermometer. Equation 


Sus rewritten as: 


Z jez 


u ER 2 

DL OmU I) S 
Toy is the parameter obtained by measuring t from the son c 
thermometer. Being a virtual temperature, this measure- 


ment is subject to humidity effects. 


B Discuss tonto Possible Errors in Wind Velocity 
and Me ineraturemleasurements 


Errors are introduced by humidity, temperature and 
atmospheric pressure which require correcting the readings 
from the sonic anemometer/thermometer. Detailed discussion 
of errors resulting from misalignment of the array, line- 
averaging,finite spatial resolution of the array and noise 
in the detector has been described by Kaimal et al (1971). 
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Me Ms uac on (3) and equation (10), the wind 
velocity measured will be a function of temperature, humid- 
ity and atmospheric pressure. From eguation (12), however, 
the measured temperature will depend on humidity and atmos- 
pheric pressure only. For hot humid days tne effect on 
NENE y cannot be ignored. Humidity significantly effects 
the absolute temperature but is negligible as far as fluctua- 
tions in temperature are concerned provided the fluctuations 
in humidity are small. 


The ratio K, between the wind velocity V (at the 








20 
/ 
temperature Too + 293K or 20C and e/y = 0)” and the wind 
velocity Vers y (at the temperature Toy) is 
_ Vasv ` Toy DONNE Du 
en oa see 
20 20 

where Boy = sound virtual temperature measured by the sonic 
thermometer in degrees C. Any wind velocity or temperature 


measured by the sonic anemometer/thermometer in an atmos- 
phere other than 20C and zero percent relative humidity 


Ee Corrected by a factor K (Figure 11A, 11B). 


2 2 
The error introduced by the assumption of C >> V 
is 
2 
iv = — (14) 
C 


and is such that for wind speeds of less than 10 m/sec it 


is zero (Figure 12). 


Sonic anemometer/thermometer circuitry is set for a base 
temperature of 20C and zero percent relative humidity. 
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Figure 11A. Temperature Compensating Constant. 
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Figure 11B. Humidity Compensation Graph. 
Kaijo Denki (1969) 
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Lc, aan Sequence 


A more detailed 
analysis by Soumie and 
Businger (1958) showed 


that velocity and tempera- 





mire within 0 10 20 30 40 50 
V (sec) 


certain limits, were pro- 


Portsonal to the fluctua- Figure 12. Error by Wind 
velocity Kaijo Denki (1969) 
tions in the difference 


and the sum of the time in- 

tervals respectively. A discussion of the method of obtain- 
ing these time intervals is of use in order to understand 

the setting up and calibration of the sonic anemometer/ther- 
mometer. Only the vertical transmitter-receiver array (that 
used to measure the vertical velocity component and tempera- 
ture) needs to be considered in order to understand the 
triggering and pulse measuring system used in the instrument. 


Pulses are emitted simultaneously from TR and TR.» 


at time T, and reach their respective receivers R and R 


Wl 2 


attes stimme intervals of ft and t. seconds (Figure 13). 


wl 2 


The time difference At (ENE - t is converted 
W w2 wl 


by a wind analog card (within the w unit) and changed into 
a voltage representing the vertical component of the wind 


velocity The output of Ro: drives the delay card Di which 


generates the delay time t This output synchronizes the 


pu 


generation of the next pulse or oscillation at time T' which 
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teach their receivers after time intervals SN and Ee 


seconds. Therefore, the pulse spacing is tyl and tni and 
the total time from the emitting of the first pulse, T, to 


is t + t 


the receiving of the second pulse, T', by R > wl Di 


ai 
S w2 


Iinsetermssorsthe speed of sound in still air, C, and 


the separation, %, of the transmitter and receiver, 


22 
t e 
E D DA vog w2 C S E (15) 


As previously discussed, the temperature measurement is 
sampled at half the rate of the wind measurement (220 times/ 


sec as compared with 440 times/sec). Therefore, the time, 


t= t ÓN is related to temperature. Since At, changes 
wl w2 D 

very slowly with temperature, a time tna generated by 

another time delay circuit is subtracted from t + Loi: The 


time difference, At, = t + t 


9 DI ^ tp2» is then supplied to a 


temperature analog card which changes it into a At, (the 
mean temperature found on the 0 dial of the temperature 
unit. The dial is operated manually). The time difference 


is also changed into At' the fluctuating component, which 


6? 
is changed into a voltage that appears on the O ov meter 
(Figure 14). 


The time delay t is adjustable (longer at lower 


D2 
temperatures) since it is preferable to keep Atp small. If 
At, is made too small, however, a sluggish limiting of the 


temperature signal will occur--the effect showing as a 


rounding to flattening of the high temperature peaks. 


SRI 





" 
DL — Bëss 5 


Time 





Figure 13. Timing Sequence Schematic. Kaijo Denki (1969) 
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Figure 14. Temperature and Time Relationship. 
Kaijo Denki (1969) 
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III. DATA ANALYSIS 


A. GENERAL 

Turbulent processes appear as erratic changes in tem- 
perature and in direction and speed of the wind. A know- 
ledge of the nature of turbulent processes in the atmosphere 
(fluctuations of such parameters as wind velocity, tempera- 
ture and moisture for example) can lead to a better under- 
standing of heat, energy and momentum transfers which occur 
over the ocean. 

The mechanics of turbulence are so complex that it is 
almost traditional now to use statistical methods through 
which an understanding of the mean characteristics can be 
obtained if not in detail. One statistical method which 
has been used with some considerable success is spectral 
analysis. From spectral analysis the distribution of energy 
of a fluctuating, turbulent process (e.g. wind velocity, 
temperature) can Ne plotted as a function of frequency. 

Analog measurements were recorded of wind velocity and 
temperature fluctuations from the sonic anemometer/thermo- 
meter,hot wire anemometer and a one mil thermocouple then 
digitized and then spectral analyzed by a Fast Fourier 
Transform method. A general discussion of procedures used 
in the analog-to-digital conversion and spectral analysis 
will be given in the next sections. More detailed discus- 
Elo niko (Mr rocdedures and theory appear in McKendrick (1971) 


ame zer and Wilson et al (1969). 


— =— 
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De KOLMOGOROV'S HYPOTHESIS 
l. General 

An important consideration in turbulence theory is 
whether the turbulence is anisotropic or isotropic. Aniso- 
tropy exists for large scale eddies and isotropy is thought 
to exist for small scale eddies. Isotropy was described by 
Lumley and Panofsky (1964) as follows: 

Mansoor op Ler eld is one in which there is no 
preferred direction. Thus all average functions de- 
scribing the statistics of such a field must remain 
unchanged regardless of a rotation (or reflection) of 
coordinate system. For example, the spectrum of a 
function, (k), homogeneous in three spatial directions, 
which is also isotropic, cannot be a function of the 
separate components of the wave number vector, for a 
rotation to a new coordinate system would change those. 
Man be ashunetcion only of the length of the vector, 
Noe cias the only quantity characterizing 'k' which 
does remain unchanged under a rotation." 

Kolmogorov postulated that at high wave numbers 
turbulence should be isotropic even though it is generated 
from and embedded in anisotropic eddies of low wave number.. 
ie the Reynolds number (ratio of inertial forces to viscou 
forces in the flow) is high enough, turbulence will adjust 
until a statistically steady state, independent of initial 
Ge nuditi von psa Li tai ned WWfurbulence is adjusted through 
inertial transfer (forces due to fluctuating velocity gradi- 
ents) and viscous dissipation. This steady state consists 
of energy transfer by inertial forces into the high wave 


number range that is egual to energy lost through viscous 


dissipation. This predicted behavior of velocity spectra 


25 





oe Ma m Mize dv nr inter pretaUlon and comparison of the sonic 
anemomemeter, hot wire anemometer and one mil thermocouple 
data. Predietion by Kolmogorov's hypotheses are discussed 
in more detail in the following paragraphs. 
2. Kolmogorov's First Hypothesis: Equilibrium Range 

At high Reynolds numbers "small scale" eddies are 
MOCT opic and in universal statistical equilibrium. Their 
properties depend on dissipation (€) and on kinematic 
eco Si ty MV) and not directly on time. Motion associated 
with the equilibrium range is uniquely determined by € 
and vV, 

Kolmogorov concluded that the energy spectrum 


[E(k)], in units of energy per unit wave number, is 


E (k) k 
Cer: SE GS 
(ev?, 1/4 Ko 
where f (E- = undesignated function 
S 
a _ (27) 1/4 
ko € 
k z characteristic wave number 
S 


3.  Kolomogorov's Second Hypothesis: Inertial Subrange 


If the Reynolds number is large enough, the energy 
containing eddy wave number region and the "dissipation" 
wave number region should be well separated in wave number 
Space. There should exist a range of wave numbers (eddy 
sizes) on the low "k" side of the dissipation range (but at 
higher wave number than energy containing wave numbers) in 


which the universal statistical steady state holds 
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[$ E(k) = 0] but where there is negligible dissipation of 
Bine mo MU un el from Lumley and Panofsky (1964)]. This 
range of wave numbers where the dominant feature is energy 
fix by inertial forces is called Sie inertial subrange. 
Since viscosity is not important in this wave number 
range, statistical properties are determined only by dissi- 


pation rate (€) and wave number (k). Kolmogorov stated that 


energy spectrum should have the following functional rela- 


tion 
E(k) = C EI „>43 
€ = dissipation of turbulent 
kinetic energy 
k = wave number 
C = constant of proportionality 


Since frequency is directly proportional to wave 
number, a log-log plot (base 10) of energy spectra versus 
frequency would yield a curve with a -5/3 slope if the spec- 
tral estimates are within the frequency band of the inertial 


subrange. 


C. ANALOG TO DIGITAL PROCEDURES 
1. General 
Analog data recorded in the field were digitized 
through use of a Hybrid computer en The Hybrid systen 
consists of an analog computer, COMCOR Ci-5000, electron- 


ically connected to a digital computer, XDS-9300. A program 





This system is located in Spanagel Hall, Naval Postgraduate 
School, Monterey, California. 
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puc. loredsaestusineo. using the XDS-9300 and the Ci-5000 
analog computer. An analog patchboard and a logic board, 
which are plugged into the Ci-5000 were also available. 
Digitization was performed in the octal number base onto 
seven track magnetic tapes. 

2. Considerations Prior to Digitization 

Considerations taken into account prior to digiti- 
zation were noise, D.C. level, analog signal strength and 
Portal analog data record length, T (corresponding to time). 
All of these were functions of collection technique. 

Noise that contaminates the desired signal is 
usually of high frequency and normally can be low pass fil- 
tered during the data recording procedures. This is the 
Preferred method. However, filtering can be done in the 
Bes computer facility just prior to digitizing. In this 
study, the analog data were filtered at the Hybrid computer. 
Data from the four sonic anemometer/thermometer channels 
were filtered at 31.25 Hz and the hot wire anemometer data 
and one mil thermocouple data were filtered at 100.0 Hz. 
(Filter frequencies should equal one-half of the digital 
sampling rate.) 

The D.C. level was removed from ali data during 
EIN Li bE were Necessary that the D.C. level be 
removed from the data during data recording, a band pass 


filter with a low frequency cut off of less than one hertz 


could have been used. This would preserve some of the low 
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Peeduenretes whilewedliminatine the high frequency noise and 
the D.C. level. 

All analog signals were amplified in order to take 
full advantage of the range of the Hybrid computer system. 
Ke erd systenerowable to dipitize signals of + 100 volts. 
Channels A and B of the sonic anemometer/thermometer and 
the hot wire anemometer were amplified by 50; Channels W 
and T of the sonic anemometer/thermometer and the thermo- 
Couple were amplified by 100. Amplification factors were 
Haken ito account, in this study, in the Fast Fõurier Trans- 
form program. 

Consideration had to be given to obtaining an opti- 
mum analog data record length, T, during the field experi- 
dents InWopntimizing this total record length for the 
frequencies of interest, the relationship between T and tlie 
Ella pe of thesconvolved wave forms (Fourier Transform) must 
DN onsidered^ "Jones (1971)9described fully the effects of 
record lengths on the shape of the convolved wave forms. 

In general there is Sonel optinumirecord length where the 
difference between the convolved spectrum and true spectrum 
very sma lN Although Jones (1971) did not find an opti- 
mum record length; he found that increasing a record from 
140 seconds to 360 seconds had no significant effect on 
spectral shape between 10 Hz and 1000 Hz. Data collected in 
diis study had a record length of 30 minutes. 

Hoeggen Mresolutlon is also tied directly to 


the length of the analog record. According to Rayleigh's 
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eo. Sad to destinguish low frequencies (f) in 


the interval 0 < f < 1/2T (where T < analog record length 


in seconds). When dealing with a finite record the original 
Signal is truncated. A long record narrows the frequency 
resolution. Record lengths of 30 minutes save a low fre- 


guency resolution of approximately 0.0003 Hz although sta- 
eco reliability would not be expected until 0.003 Hz. 
Another factor which had to be considered prior to 
digitizing was that of aliasing.  Aliasing is dependent on 
the folding (Nyquist) frequency which is 1/2 the sampling 
frequency. In order to prevent the'folding back" of energy 
from higher frequencies into the lower frequencies of 
interest, a sampling rate greater than twice the highest 
freguency desired must be chosen. The folded back energy 
Will be centered around the Nyguist freguency. Sonic anemo- 
meter/thermometer data from Experiment I was run at two 
ditferent sampling rates. Sampling rates of 390.5 Hz and 
deee ave aliasing frequencies of 195.25 Hz and 31.25 Hr? 
respectively. GE Wire anemometer data and the thermo- 
couple data from Experiment II were sampled at 200.0 Hz which 
gave an aliasing frequency of 100.0 Hz. Experiment LI sonic 


anemometer/thermometer data were sampled at 62.5 Hz only. 


D. CONVERSION OF SEVEN TRACK TO NINE TRACK 
Digitization of analog data on the above system produced 
a seven track tape in the octal number case. To utilize 


this data in the IBM-360/67 computer system, digitized data 
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had to be converted to a nine track tape in the hexadecimal 
number base. This conversion process was done using an 
available program, CONVERT. Complete listings and Job Con- 
trol Language (JCL) card setups for the CONVERT program can 
be found in theses by McKendrick (1971) and Jones (1971), 


and in a report by Raney (1971). 


ER FOURIER" COEFFICIENT COMPUTATIONS 

Fourier coefficients were obtained by using the FTOR 
program, one of two UBC time series programs currently on 
file in the IBM-360/67 subroutine library at the Naval Post- 
graduate School. 

Transformation of digitized signal data from the time 
domain into the freguency domain was done using a finite 
form of the Fourier Transform. Data were read off the con- 
verted nine track tape, a record at a time. The Fourier 
coefficients were computed and were then stored on a sepa- 
mate tape. 

The longer the record, the closer together in frequency 
will be the spectral estimates because the value of the 
frequency is given by (n/T). That is,the discrete frequency 
bandwidth is inversely proportional to the record length. 

Due to the prohibitive amount of computer time used to 
compute coefficients, the direct finite Fourier Transform 
of a series of data points was not used. The Fast Fourier 
Transform of Cooley and Tukey (1965) was used and it is a 


Wise nok the direct method. The Hast. Fourier Transtorm 
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is incorporated into the FTOR program resulting in consider- 
du MossNcompurerctime for Fourier coefíicient computation. 
This method yields an improvement in excess of five orders 
of magnitude in computational time (McKendrick 1971) over 
the direct method. 

The FTOR program transforms a maximum of 217° (8192) data 
points. This maximum is set by dimension limitations in the 
program. Turbulence data normally includes many more data 
points than this so the data points are read in as records 
consisting of integral powers of two (up to a maximum of 


13 
2° ) and these records are transformed consecutively. 


F.  UBCSCOR PROGRAM 

Once the Fourier coefficients have been computed, the 
second UBC program, SCOR, converts them to power spectral 
densities. 

Spectral values E. (n) are computed from the complex 
Fourier coefficients by forming the product of the complex 
coefficient eu) = SE - i DM with a conjugate 
[A (n) - a, (n) + i b, (n)] where j and k could be the same 
or different. | 

Spectral values were obtained by the following computa- 
tions 

Specetremwdensity2 j- = k, 
* 


A. (n)*A, (n) 


= Ek e = 18 
Een) py (a Sg bo ) 2 


where An = 


aja c 
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Cross spectrum: j % k, 
* 
A. (n) * A, (n) 


EEUU, = An 


li 


T 

se + — 4 = 

ma NO usc APE, 

The real and imaginary parts of the cross spectrum are 


the cospectrum and quadrature spectrum respectively and 


were denoted as; Cospectrum: 


T 
a = (ac, + bd) 7 


Quadrature spectrum: 
Da a) > (orco E aid? Då 

The spectra and cospectra are computed for each fre- 
quency in successive records and then these records are 
averaged consecutively. The program also divides the spec- 
tral densities into 32 bandwidths and averages them across 
the bandwidths as well as along the records to obtain a 
mean and a variance of the quantities. 

In the SCOR program one can specify the number of 
channels to be analyzed, the number of records in a file, 
the start point for analysis, selection of spectra and/or 
cos peca band dela desired (logarithmic or linear), the 
type of graphic plot desired and for this study, machine-set 
Ccaordinates for the graphic output. Programmer-set-co- 
ordinates also can be used. The program computes, the 
spectra and/or cospectra for each record in sequence until 


all records requested have been analyzed or the "end-of- 


file" mark is reached on the tape. 
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TERE ONE PDERIMENES 


AS PRELIMINARIES 

The sonic anemometer/thermometer used in the field ex- 
periments had been in storage for a period of two years, 

The first consideration was to thoroughiy check the sonic 
ie moneter/thermnometer electronically. It was checked out 
meae cordance with manufacturer's instructions [Kaijo Denki 
(1969),Appendix D and E which included voltage and wave form 
checks at certain critical points]. A large number of tran- 
sistors in both the power supply and junction box had to 

be replaced. 

Once the sonic anemometer/thermometer was checked elec- 
tronically, wind tunnel tests were performed. Simple, 
operational tests in a controlled environment were made to 
insure that the sonic anemometer/thermometer was measuring 
wind speed and temperature correctly. Once satisfactory wind 
tunnel EE were made, field experiments were 


conducted. 


BE UDESREOLIVESVOF EXPERTIENTS 

The objective of the field experiments was to gather com- 
parative turbulence data from the sonic anemometer/thermo- 
meter,the hot wire anemometer and a one mil thermocouple. 
All three sets of data were then similarly spectrally 


analyzed. 
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C. ESBERILMENTS 
l. Installation of Instruments 

The procedures described in this section applies to 
both experiments. All equipment was installed at the Naval 
Postgraduate School beach lab,indicated in Figure 16 and 
Figure 17. The sonic probe, hot wire probe and one mil ther- 
mocouple were mounted on a pole at a height of three meters 
(Figure 18) and positioned so that all sensors were facing 
directly into the wind flow which was blowing off the water. 
During the experiments, all probes were kept into the wind 
such that the velocity measured on channels A and B were 
equal. A hand held wind vane was used to check on the wind 
direction in order to maintain proper probe orientation. 

All recording devices,power supplies and amplifiers 
were installed inside the beach lab (Figure 19). Analog 
data were recorded on a SANGAMO FM tape reccrder model 3562 
at a speed of 1-7/8 inches per second. This recorder uses 
one inch magnetic tape and has 14 tracks available for data 
and two edge tracks available for voice recording. 

The hot wire anemometer used for velocity compari- 
sons was a TSI model 1050. A copper constantan, unsheathed, 
fine gauge Omego Engineering Company (SSTSSP-020E-6) thermo- 
couple was used for temperature measurements. 

DEMNM emi menim 14 December 1972 

Wind speed and temperature measurements were mad? 

during the late afternoon on 14 December 1972. isibility 


was unlimited. Winds were from the west to northwest, off 
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the ocean, at an average speed of 2-4 m/sec. Mean air tem- 
Eersceure tor the period (1430-1700 PST) was 52°F., 

Three separate 30 minute analog data records were 
Her or demo WO0OME o 1550, 1546 to 1616 and 1620 to 1650 PST). 
Data on all channels were recorded without amplification. 
During the measurements, spikes were observed on all channels. 
The spikes were found, after this experiment, to have been 
caused by faulty transistors in the junction box.  Examples 
of these spikes appear in the strip chart recording in 
Figure 20. 

Hot wire anemometer and thermocouple data were not 
recorded because the record amplifier on the hot wire channel 
of the tape recorder failed and the thermocouple malfunc- 
tioned. 

The three separate analog data records from the 
sonic anemometer/thermometer data were digitized and spectral 
analyzed. 

Ore nen - 27 February 1973 

This experiment was similar to experiment I except 
that hot wire anemometer data and thermocouple data were 
successfully recorded. Channel B of the sonic anemometer/ 
thermometer was not operating because of a short in the con- 
nector cable on that channel. This malfunction did not 
affect the other channels. 

Again, wind speed and temperature data were obtained 


during three separate 30 minute periods (1030 to 1100, 1345 
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pom oe nd 1530 to 1600 PST). It was an overcast sky, 
Visibility was approximately six miles. Winds were from the 
west, off the ocean, with an average speed of 4-8 m/sec. 
Mean air temperature ranged from 58°F in the late morning 
period to 60°F in the early afternoon period. 

During the recording periods, all but one of the 
data channels were amplified using Preston 8300-XWB ampli- 
fiers. Ihe amplification factors were two for the sonic 
W and T channels, one for the sonic A channel, twenty for 
the hot wire channel and one thousand for the thermocouple 
channel. The hot wire data were filtered before recording 
using a Krohn-Hite model 3750 bandpass filter low set at 
0.2 Hz and high set at 200.0 Hz. The three separate 30 
minute analog data records were digitized and spectrally 


analyzed. 
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V. RESULTS AND CONCLUSIONS 


A. GENERAL PERFORMANCE OF SONIC ANEMOMETER/THERMOMETER 

In general the sonic anemometer/thermometer satisfac- 
torily measured the wind and temperature fluctuations during 
the field experiments, but it did not operate at its peak 
potential because of various electronic problems. For 
example, spiking re òn allchannels during both ex- 
periments I and II and a short occurred in channel B during 
experiment II. These could have been a result of the instru- 
ment being idle for two years. 

Spiking proved to be a continuous problem during the 
field experiments even though numerous transistors in the 


system had been replaced. 


D. SPECIRA PLOIS 
l. General 

Sonic anemometer/thermometer data recorded on 14 
December 1972 (experiment I) were digitized at a sampling 
palio Mo i 90% 5 Hz (aliasing frequency 192.25 Hz) and also 
IO eo aliasing frequency 31.25 Hz. 

No editing or filtering was performed on any of the 
channels during digitization at the higher sampling rate. 
The higher sampling rate (390.5 Hz) enabled the examination 
of a frequency band from © Hz tO approximately 100 Hz. 

The analog data tape was edited prior to digitizing 


at the lower sampling rate. This led to the elimination of 
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A nep ento which had the heaviest spiking. The 
remaining two time periods were digitized at 62.5 Hz and 
filtered at 31.25 Hz. The lower sampling rate shifted the 
frequency band examined (approximately 0.5 Hz to 40 Hz) 
towards the lower end of the frequency spectrum. The spec- 
tral results indicate both the energy containing frequency 
band as well as the lower frequency end of the inertial sub- 
range. 

Hot wire data were not taken during experiment I 
because of a recording amplifier failure. Two samples of 
hot wire data taken on a day that had similar synoptic 
characteristics (air temperature and wind) as the day of 
experiment I are included for comparison purposes. They 
are considered representative samples of data which would 
have been taken by a hot wire anemometer system on that day. 
These results represent a frequency band of from 0.01 Hz 
to approximately 100 Hz. The data for both of these results 
were digitized at a rate of 200.0 Hz (aliasing frequency 
100.0 Hz) and Hi Peed ac MODNO H z% 

On 28 February 1973 (experiment II), sonic anemo- 
meter/thermometer, hot wire anemometer and one mil thermo- 
couple data were collected during the three 30 minute time 
periods, The sonic data were digitized at a sampling rate 
emo ome onde: titeredq at 32.25 Hz. The hot wire and 
thermocouple data were digitized at a sampling rate of 
200.0 Hz and filtered at 100.0 Hz. Sonic channel B was 


inoperative. 
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2. Sonic Anemometer/Thermometer Data Sampled at 


oe o December 1972 (Experiment I) 


Figure 21 to Figure 32 are for the three 30 minute 
periods from 1500 to 1700. A -5/3 slope exists on the log- 
ori ao cia the frequency band from 1.0 Hz to 10.0 Hz for 
TE hannels. Spectra from channels A, B, W have -5/3 
slopes which continue beyond 10.0 Hz but the T spectra 


Ses drop off more rapidly than -5/3 before 10.0 Hz. 


S. Sonic Anemometer/Thermometer Data Sampled at 
DNI l4 December 1972 (Experiment I) 


Figure 33 to Figure 40 are from the two edited 
periods of data which were filtered and then sampled at the 
lower rate, Spectra from channels A and B in both periods 
Mente —5/3 slopes within the 1.0 Hz to 10.0 Hz band 
with some departure around 8.0 Hz. For both periods, W 
channel spectra have -5/3 slopes up to 6.0 Hz with rapid 
deterioration beyond. For the second period, the T spectrum 
Da drop off from a —5/3 slope at approximately 4.0 Hz but 
moi ne third period, the T spectra has a definite -5/3 


slope throughout the 1.0 Hz to 10.0 Hz frequency band. 


4. Hot Wire Anemometer Data Sampled at 200.0 Hz 
TALA 197/3 


Figure 41 and Figure 42 are from hot wire anemometer 
data taken on 10 January 1973. Two periods are discussed 
and each period was approximately 50 minutes in length. A 
distinct -5/3 slope is found in the frequency band from 


IO ap pro: imately" 30.0 Hz for both periods of data. 
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Figure 21. Sonic "A" Channel Data Observed 14 December 
1972 lst Period (Experiment I) High Sample Rate. 


57 





€ HERTZ) 


PISO 





LOGLO FREQUENCY 


Figure 22. Sonic "B" Channel Data Observed 14 December 
1972 lst Period (Experiment I) High Sample Rate. 
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Sonic Channel Data Observed 14 December 
Ie mod (Experiment I) High Sample Rate. 
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Figure 25. Sonic "A" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 26. Sonic "B" Channel Data Observed 14 December 
72 nd Period (Experiment I) High Sample Rate. 
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Figure 27. Sonic "W" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 31. Sonic "W" Channel Data Observed 14 December 
1972 3rd Period (Experiment I) High Sample Rate. 


67 





£ HERTZ) 


> A 
=e 8 






(TABU 


r- 
— 
LI 


—a 9 


1.0uz 
(å 
E 
1 Pr 

3 E 

Sé NU 
I ^ | 100uz 
| | E 


ZU 3 


Sle se 


2 





Es 
p 
J3 
= 
m 
FE 
C) 
ca 


1 


| 
mole 


Figure 32. Sonic "T" Channel Data Observed 14 December 
1972 3rd Period (Experiment I) High Sample Rate. 
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Her ene VW Sonic "TI" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) Low Sample Rate. 
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Figure 41. Hot Wire Data Observed 10 January 1973 
Ist Period. 
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Figure 42. Hot Wire Data Observed 10 January 1973 
Dude Perrod. 
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5. Sonic Anemometer/Thermometer Data Sampled at 
Ono Hz 725 February 1973 (Experiment II) 


Figure 43 to Figure 5l are spectra from the sonic 
anemometer/thermometer data recorded on 28 February 1973. 
Channel B was inoperative during the period.  Spectra from 
channels A and W — a definite -5/3 slope in the 1.0 Hz 
to 10.0 Hz frequency band. A -5/3 slope in a very narrow 
band of frequencies of approximately 1.0 Hz to 4.0 Hz 


appears in the T spectra. 


6. Hot Wire Anemometer Data Sampled at 200.0 Hz 
- 28 February 1973 (Experiment II) 


Figure 52 and Figure 53 are hot wire spectra from 
the second and third periods recorded on 28 February 1973. 
The first period of data was lost due to a faulty magnetic 
tape. Spectra from both periods have definite -5/3 slopes 


faethe region 1.0 Hz to approximately 91.0 Hz. 


7. Thermocouple Data Sampled at 200.0 Hz - 
Me busty 19/3 (Experiment Il) 


Figure 54 and Figure 55 are thermocouple spectra 
from the second and third periods recorded on 28 February 
1973. The first period of data in this case was also lost 
due to the faulty magnetic tape. From the second period 
plot the -5/3 slope starts at approximatelv 4.4 Hz and goes 
ucc ug) 0sHz.— From the third period plot, a 


definite -5/3 slope exists in the range 1.0 Hz to 10.0 Hz. 


Gi CONCLUSIONS 
The sonic anemometer/thermometer system proved to be a 


reliable, rugged and easily installed instrument for the 
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Figure 43. Sonic "A" Channel Data Observed 28 February 
1973 lst Period (Experiment II). 
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Figure 44. Sonic "W" Channel Data Observed 28 February 
1973 lst Period (Experiment II). 
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Figure 47. Sonic "W" Channel Data Observed 28 February 
Sonda Period (Experiment II). 
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Figure 48. Sonic "T'" Channel Data Observed 28 February 
1973 2nd Period (Experiment II). 
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Eure 749, Sonic A’ Channel Data Observed 28 February 
re ra Period (Experiment II). 
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Figure 51. Sonic "T" Channel Data Observed 28 February 
arg period (Experiment II). 
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Figure 52. Hot Wire Data Observed 28 February 1973 
2nd Period (Experiment II). 
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Figure 53. Hot Wire Data Observed 28 February 1973 
Ena perjod (Experiment II). 
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Eigure 54,  Thermocouple Data Observed 28 February 1973 
2nd Period (Experiment II). 
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Figure 55.  Thermocouple Data Observed 28 February 1973 
Or dW er iod (Experiment II), 
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measurement of the three components of the wind and also air 
temperature. Its freguency response was found to be from 
ISE" to approximately 12.0 Hz. 

The sonic anemometer/thermometer measurements did 
not appear to be affected by weather elements and because it 
uses the speed of sound in air it has inherent calibration. 
The response of the sonic anemometer/thermometer is linear 
and so is its sensitivity. It is thus able to be used to 
describe vertical and horizontal transfers of energy, heat 
and momentum. 

The -5/3 slope of Kolmogorov's inertíal subrange 
was easily recognizable on spectral plots of turbulence data 
measured by the sonic anemometer/thermometer. Thus the com- 
putation of the rate of dissipation of kinetic energy con- 
stant can be made which in turn means that heat and momentum 


transfer computations can be made. 
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APPENDIX A 


Definition of Terms Used in Analog-Digital Data Processing 


The terms used in the analog-to-digital process are de- 
fined here based on current usage in the Hybrid computer 
center and based on that utilized by the three time serles 
programs for the spectral analysis of digitized data. There 
are some contradictions of terms between McKendrick (1971) 
and Jones (1971) and between what is currently used. 

The word "record" is often used in both analog and digi- 
tal discussions as meaning a finite length of data in 
seconds. In the analog sense, "record" means the entire 
imo VO data which is to be digitized., In the digital 
sense, "record" means a short length of data based on the 
sampling rate and the number of samples per record selected 
by the programmer. In McKendrick (1971) and Jones (1971), 
digital record" is referred to as "block". The word 
"block" is not used in this thesis as "record" is what is 
Meco alll three programs (although the print outs use 
blocks). A digital "record" is limited in the number of 
samples it can contain because of Hybrid computer memory 
limitations., This limit is that the number of samples times 
the number of channels simultaneously digitized must be 
egual to or less than 6144 channel-samples. As an example, 


record 
consider the following hypothetical problem: 
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Analog record length = 3,000 secs 
Number of analog channels containing data = 5 channels 


Sampling rate = 100 samples 
sec 


Number of channels digitized 


simultaneously = 5 channels 
Samples per record selected = 1024 samples 
record 


1024 samples x 5 channels = 5120 channel-sample 


record — (601044 
record E 
samples l sec n seconds 
Dr record 100 samples oboda record 


Digital record length < 10.24 seconds 
record 


Note that the value for samples for record was chosen as 
an integral power of two , i.e. E The reason is that the 
maximum number of data points (samples) that can be stored 
in the IBM 360/67 computer memory under the FTOR program are 
8192 or nos samples. Therefore, for efficiency purposes tae 
FTOR program was designed to read in samples per record wt .ch 
contain integral powers of two. 

Omemorhier commonly wsed word is. fille." A "file" is a 
finteestenegth of digital data containing any number of 
records from one on up. The length of a "file" is predeter- 
ue oda ne the digitization phase and is based on the 
number of records desired per "file." An "end of file" 


mark is placed on the seven track tape after the desired 


number of records have been computed. 
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